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Heavy flavour production is one of the key components of the HERA II physics programme. While most of the 
results presented use leptons or the reconstruction of charmed mesons to identify heavy flavour production, both 
the HI and ZEUS experiments now have working microvertex detectors that are being used more and more. In 
this talk I will summarise a selection of the recent results obtained by the two collaborations. 



1. Introduction 

The study of heavy flavour production at 
HERA provides an important test of pertur- 
bative Quantum Chromodynamics (pQCD) and 
also valuable information for the measurements to 
be made at the LHC. I will discuss a small selec- 
tion of the many measurements of heavy flavour 
production that have been made at HERA, con- 
centrating on the more recent results. 

HERA running started in 1992 and the ac- 
celerator stopped running at the end of June 
2007. HERA collided 27.5 GeV electrons or 
positron^] with 920 GeV (820 GeV until the end 
of 1997) protons, giving a centre-of-mass energy 
of 318 GeV. A long shutdown in 2000 and 2001 
was used to upgrade the machine and the detec- 
tors, with the aim of increasing the luminosity by 
about a factor of 4. The period up to 2000 is usu- 
ally called HERA I and after 2001 HERA II. By 
the end of the running both of the colliding beam 
experiments, HI and ZEUS, had collected about 
0.5 fb- 1 of data. 

Several kinematic variables are used to charac- 
terise the ep scattering process: 

• Q 2 , the negative squared four-momentum 
exchanged at the electron or positron ver- 
tex; 

• x, the Bjorken scaling variable; 

• y, the inelasticity; 

• W, the invariant mass of the hadronic final 
state. 

1 Hereafter unless explicitly stated both electrons and 
positrons are referred to as electrons. 



The measurements are usually separated into 
the deep inelastic scattering (DIS), Q 2 > 1 GeV 2 , 
or photoproduction, Q 2 < lGeV 2 , regimes, de- 
pending on whether the scattered electron is de- 
tected in the in the main calorimeter or not. 

The main production mechanism for heavy 
quarks is the so-called Boson Gluon Fusion 
(BGF) process which is illustrated in Figure [I] 
In practice higher order contributions also have to 




Figure 1. Feynman diagram of the boson gluon 
fusion process. 



be taken into account. The photon can fluctuate 
into a qq pair and one of the quarks then partici- 
pates in the hard interaction (resolved photopro- 
duction); or in so-called excitation the scattering 
from charm or beauty can take place with intrin- 
sic charm or beauty inside the proton or photon. 
Monte Carlo models usually include these pro- 
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cess in addition to the direct boson-gluon fusion. 
These higher order processes are also sometimes 
referred to collectively as non-direct. 

Next-to-leading oder (NLO) QCD calculations 
exist for heavy quark production at HERA. 
These are implemented for photoproduction in 
the FMNR programme [IE] and for DIS in the 
HVQDIS programme [3 j. The programmes in- 
clude simple independent fragmentation of the 
b or c quark, but on their own are not able to 
give predictions for correlations between final- 
state particles. The FMNR programme has been 
interfaced [4] to the Pythia Monte Carlo and its 
predictions have been used for studies of dimuon 
final states. 

Most heavy flavour measurements rely on the 
central tracking detectors and are helped signifi- 
cantly by the presence of a microvertex detector. 
HI had such a detector for some of the HERA I 
running, while both ZEUS and HI had such de- 
tectors for the HERA II running period. 

Several different methods have been used by 
the collaborations to identify the production of 
heavy quarks. Each of them has its advantages 
and disadvantages and often cover different kine- 
matic ranges. Traditionally the identification of 
D* mesons and semileptonic decays have been 
used most often. With the advent of microver- 
tex detectors lifetime information is being used 
more and more often. This works best for events 
with high energy jets, while tagging with semilep- 
tonic decays to electrons or double tags enables 
one to go to lower transverse momenta. 
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2. D* Production 

The HI collaboration has made a series of mea- 
surements of D* production both in the photo- 
production [5] (96 pb -1 of data taken in 2006 
and 2007) and the DIS [6] (347 pb" 1 from the 
HERA II running period) regimes. For these mea- 
surements they make use of their new Fast Track 
Trigger, which enables events with D* candidates 
to be selected early in the trigger chain. In both 
photoproduction and DIS clear signals are seen 
as illustrated in Figure [2j The kinematic cuts are 
indicated in the figures. 

Differential cross-sections as a function of a 



Figure 2. AM = m(K7nr) —m(K7r) distributions 
for the photoproduction (top) and DIS (bottom) 
samples. The points show the data, the curves 
in the lower figure, the results of a fit to the dis- 
tribution. The dotted line shows the background 
shape. The kinematic cuts are indicated. 
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wide range of variables have been determined. 
The cross-section as a function of Q 2 is shown 
in Figure [3] The measurements are compared to 
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Figure 3. D* production cross-section as a func- 
tion of Q 2 . The data (points) are compared to 
NLO QCD predictions with two different PDFs. 
The lower plot shows the NLO QCD prediction 
scaled by the ratio of the data to NLO QCD vis- 
ible cross-sections. The range of variation of the 
quark mass and factorisation and renormalisation 
scales are indicated in the figure. 



the predictions of the HVQDIS programme us- 
ing two different Parton Distribution Functions 
(PDF). The different PDFs show a very similar 
behaviour as a function of Q 2 . 

In contrast the cross-section as a function of the 
pseudorapidity, 77, shows a clear dependence (see 
Figure [4|, demonstrating the sensitivity of charm 
production to the gluon structure function. 

Making the same comparison for photoproduc- 
tion one sees that the NLO QCD prediction un- 
derestimates the data in the forward direction 
(see Figure |5|. Comparing the data to Monte 
Carlo predictions, the Cascade generator, which 
is based upon /^-factorisation, is able to pro- 
vide a better description of the data in DIS than 
RAPGAP with two different PDFs (see Figure^. 
However, in photoproduction Pythia, using the 
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Figure 4. D* production cross-section in DIS as 
a function of r](D*) compared to NLO QCD pre- 
dictions. For further details see the caption of 
Figure [3j 
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Figure 5. D* production cross-section in pho- 
toproduction as a function of r](D*) compared 
to the NLO QCD prediction. The band shows 
the uncertainty in the prediction. The range of 
variation of the quark mass and factorisation and 
renormalisation scales are indicated in the figure. 
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Figure 6. D* production cross-section in DIS as 
a function of r](D*). The data are compared to 
the different Monte Carlo generators as indicated 
in the figure. 



rate 6-quark events from background. The cross- 
sections as a function of the transverse momen- 
tum and pseudorapidity of the muon are shown 
in Figures [7] and [8j Both the charm and beauty 
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massless scheme for the generation of the heavy 
quarks, provides the best description of the data, 
while Cascade and Pythia in the massive mode 
both undershoot the data in the forward direc- 
tion. 

3. Beauty in Photoproduction 

The ZEUS collaboration recently reported two 
new measurements of 6-quark production in pho- 
toproduction. The first measurement uses semi- 
leptonic decays to muons [7] to identify heavy 
quark decays, while the second one uses electrons 
[8]. Dijet events are selected by requiring at least 
two jets with \rj\ < 2.5 and a transverse mo- 
mentum (muon measurement) or energy (electron 
measurement) greater than 7 GeV for the highest 
transverse energy and 6 GeV for the 2nd highest 
transverse energy jet. 

The first measurement uses 124 pb -1 of data 
collected in 2005 and requires a well identified 
muon with pt > 2.5 GeV and —1.6 < r] < 2.3. 
The microvertex detector is used to measure the 
impact parameter (S) of identified muons. This is 
combined with the transverse momentum of the 
muon with respect to the jet axis (p^ 1 ) to sepa- 



Figure 7. 6-quark production cross-section in di- 
jet photoproduction as a function of the trans- 
verse momentum of the muon. The measurement 
is compared to an earlier ZEUS publication [9 
and the NLO QCD prediction. The uncertainty 
on the prediction is indicated by the band. 



contributions are left free in the fit. The fig- 
ures also show the results of an analysis using the 
HERA I data [9] , which only used p^ 1 for separa- 
tion. For this earlier analysis the charm content 
was fixed. Good agreement between the analyses 
is seen. The predictions of the NLO QCD predic- 
tion are also in good agreement with the data. 

The second measurement identifies electrons 
from the semileptonic decays of heavy quarks us- 
ing an integrated luminosity of 120 pb -1 collected 
during the HERA I running period. Electron 
identification uses the measurement of the spe- 
cific energy loss, dE/dx, in the Central Track- 
ing Detector, CTD, the fraction of the energy 
deposited in the electromagnetic calorimeter as 
well as the ratio of the energy deposited in the 
calorimeter to the momentum measured in the 
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Figure 8. 6-quark production cross- sect ion in di- 
jet photoproduction as a function of the pseudo- 
rapidity of the muon. The measurement is com- 
pared to an earlier ZEUS publication [9] and the 
NLO QCD prediction. The uncertainty on the 
prediction is indicated by the band. 



tracking detectors. Semileptonic decays are sep- 
arated from background using p^ and the az- 
imuthal angle between the electron direction and 
the missing transverse momentum vector, Ac/). As 
illustrated in Figure [9J p^ 1 can separate 6-quark 
from c-quark decays, while Acj) separates semilep- 
tonic decays from background. 

The variables are combined using a likelihood 
ratio method, which is optimised for the iden- 
tification of electrons from semileptonic 6-quark 
decays. The distribution of the likelihood ratio is 
The distribution is fit using 



shown in Figure 10 



the expected distributions for beauty, charm and 
background to determine the fractions of events 
from each source. The fit provides a very good 
description of the data. 

The visible ep cross sections (at hadron level) 
for 6-quark and c-quark production and the sub- 
sequent semileptonic decay to an electron with 
Pt > 0.9 GeV in the range |?7 e | < 1.5 in photo- 
production events with Q 2 < 1 GeV 2 and 0.2 < 
y < 0.8 and at least two jets with Et > 7(6) GeV, 
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Figure 9. Variables used to separate semilep- 
tonic heavy quark decays from background. The 
solid line shows the distribution for electrons from 
semileptonic 6-quark decays, the dashed line for c- 
quark decays and the dotted line the background 
(Bkg). 
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Figure 10. The distribution of the likelihood ra- 
tio for electron candidates, Ac an d, in data com- 
pared to the Monte Carlo expectation after the 
scaling the predictions to best fit the data. The 
arrow indicates the region included in the fit 
(—2 InT < 10). The shaded areas show the fitted 
contributions from b quarks, c quarks and back- 
ground as denoted in the figure. 
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\rj\ < 2.5 were determined separately for y/s = 
300 GeV and y/s = 318 GeV. For the complete 
data set (96-00) the cross-sections evaluated at 
y/s = 318 GeV are 

al is = (125 ± ll(stat.)tn(syst.)) pb, 
a vis = ^ 278 ± 33( sta t.)t^(syst.)) P b - 

These cross-sections are in agreement with the 
corresponding NLO QCD predictions: 



_NLO 



NLO = (3 8 0+48) pb . 



The cross-sections as a function of the trans- 
verse momentum and pseudorapidity of the elec- 
tron are shown in Figure [TT] Good agreement 
with the NLO QCD prediction is seen. The data 
also agree well with the Pythia prediction scaled 
by a factor of 1.75 for 6-quark production and 1.28 
for c-quark production. 

4. Beauty Correlations 

The identification of both heavy-quark decays 
in an event has several advantages: the back- 
ground is reduced substantially and the kinematic 
range accessible is larger. The disadvantage is 
a significant reduction of statistics. If both b- 
quark jets can be identified, dijet correlations 
can be directly measured which probe next-to- 
leading order effects. The ZEUS collaboration 
used the HERA I data sample to select events 
with Et > 8 GeV and two identified muons [10] . 
Separating the samples according to the relative 
charge of the muons as well as their invariant 
mass allows much of the background to be evalu- 
ated directly from the data. 

The extracted cross-section is shown in Fig- 
It is compared with the NLO QCD predic- 
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tion as well as with a leading order Monte Carlo. 
Given the limited statistics it is not possible to 
say whether the NLO calculation provides a bet- 
ter description of the data. 



F| 6 and Ff 



The HI collaboration have used the impact pa- 
rameter significance to determine F$ b and F$ c for 
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Figure 11. Differential cross sections as a func- 
tion of a), c) the transverse momentum and b), 
d) the pseudorapidity of the electrons. Plots a) 
and b) are for 6-quark production while c) and d) 
are for c-quark production. The measurements 
are shown as points. The inner error bar shows 
the statistical uncertainty and the outer error bar 
shows the statistical and systematic uncertainties 
added in quadrature. The solid line shows the 
NLO QCD prediction after hadronisation correc- 
tions, with the theoretical uncertainties indicated 
by the band; the dashed line shows the scaled 
prediction from Pythia. 
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Figure 12. The dimuon cross-section as a func- 
tion of the azimuthal angle between the muons 
in dijet events. The data (points) are compared 
to the NLO QCD prediction calculated using the 
FMNRxPYTHIA interface as well as to the RAP- 
GAP prediction scaled by a factor of 1.84. The 
inner error bar shows the statistical uncertainty 
and the outer error bar shows the statistical and 
systematic uncertainties added in quadrature. 



jg10 7 
o 

,6 



Measurement of F^ c and F^ b 
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Figure 13. Distribution of the signed impact pa- 
rameter for the data (points) and the contribu- 
tions from 6, c and light quarks, evaluated using 
Monte Carlo events. 



12GeV 2 < Q 2 < 650 GeV 2 and 0.0002 < x < 
0.032 using 56 pb" 1 of data taken in 2006 [II]. 
The distribution of the impact parameter in the 
transverse plane is shown in Figure [13] 

In order to separate beauty and charm from 
background the two tracks with the most signifi- 
cant impact parameters (Si and £2 respectively) 
are selected, rejecting events where the signs of 
the impact parameters differ. The significance is 
defined as 5/a(5), where <j(5) is the error on S. 
Events with one good track are used to make the 

51 distribution; all other events are used for the 

52 distribution. 

The contents of the negative significance bins 
are subtracted from the corresponding positive 
significance bins. This yields the distribution of 
S2 shown in Figure [l4j The charm contribution 
dominates the distribution. At high significance 
the beauty contribution becomes larger. The data 
are split into bins in Q 2 and x and the contribu- 
tions of beauty and charm are determined sepa- 
rately in each bin using a least squares simultane- 
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Figure 14. Distribution of the second highest sig- 
nificance and the contributions from 6, c and light 
quarks, evaluated using Monte Carlo events. 
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ous fit to the Si and 62 distributions. The overall 
normalisation is determined by also including in 
the fit the total number of inclusive events with- 
out any cut on the impact parameter significance. 
The results of the fit in each x — Q 2 bin are con- 
verted to a "reduced cross-section" using 



a cc >,Q 2 ) 



xQ 4 



dxdQ 2 27ra 2 (l + (l-2/) 2 ) 



The reduced cross section for cc production as a 
function of x in different Q 2 bins is shown in Fig- 
ure Measurements made with the HERA I 
data in the same kinematic range are also shown 
in the figure. The results are in good agree- 
ment with each other and can be combined, tak- 
ing into account the correlated systematic errors. 
The prediction using the Variable Flavour Num- 
ber Scheme (VFNS) in the MRST04 PDF agrees 
well with the data. 

From the reduced cross-section F^ c can be ex- 
tracted: 

a c5 (x,Q 2 ) = FZ 



(l + (l-y)2) 



The measurements of Fn 



and F| 6 



The correction due to the longitudinal structure 
function, F£ c , is small. The same formulae with c 
replaced by b can be used to extract a bb and F| 6 . 

are shown in 
Figures [16] and [TTJ respectively. The rapid in- 
crease in F£ c as a function of Q 2 at low x is 
clearly visible. The F| 6 measurements are con- 
sistent with the same trend, but are less precise. 



6. Conclusions 

The many HERA measurements of beauty pro- 
duction in photoproduction are compared in Fig- 
The measurements presented here agree 
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well with the previous values, giving a consistent 
picture of 6-quark production in ep collisions in 
the photoproduction regime, and are well repro- 
duced by the NLO QCD calculations. For all the 
measurements leading order Monte Carlo predic- 
tions also describe well the shapes of the distri- 
butions. 

There is a tendency for the charm measure- 
ments to overshoot the predictions in the forward 
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Figure 15. The reduced cross- section for cc pro- 
duction in different Q 2 bins as a function of x. 
The inner errors bars show the statistical uncer- 
tainty, the outer error bars represent the statisti- 
cal and systematic uncertainties added in quadra- 
ture. The measurements are compared with those 
from HERA I and the averaged data are also 
shown. The measurements are compared to the 
MRST04 prediction. 
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Figure 16. F% c as a function of Q 2 for different x 
ranges. Also shown are measurements from the 
HI and ZEUS collaborations using D* mesons to 
identify the charm quarks. The inner errors bars 
show the statistical uncertainty, the outer error 
bars represent the statistical and systematic un- 
certainties added in quadrature. 



Figure 17. F% b as a function of Q 2 for different 
x ranges. Also shown is a preliminary measure- 
ment from the ZEUS collaboration using the data 
from 2004. The inner errors bars show the sta- 
tistical uncertainty, the outer error bars represent 
the statistical and systematic uncertainties added 
in quadrature. 
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Figure 18. Differential cross section for 6-quark 
production as a function of transverse momen- 
tum, pjn, compared to the results of previous 
ZEUS measurements as indicated in the figure. 
The measurements are shown as points. The in- 
ner error bar shows the statistical uncertainty and 
the outer error bar shows the statistical and sys- 
tematic uncertainties added in quadrature. The 
solid line shows the NLO QCD prediction from 
the FMNR program with the theoretical uncer- 
tainty shown as the shaded band. 



direction. Comparing charm quark predictions 
using different PDFs shows that the cross-section 
is, as expected, sensitive to the gluon PDF. 

In the near future a number of final HERA I 
measurements will be published. With the 
HERA II data the kinematic range of the mea- 
surements can be extended and a combination 
of different tagging methods should increase the 
precision of the measurements. The improved 
HERA II forward tracking will allow much im- 
proved studies of heavy quark production in the 
forward direction. 
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